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■"sj- ' We study the equation of motion appropriate to an inspiralling binary star 

system whose constituent stars have strong internal gravity. We use the post- 
Newtonian approximation with the strong field point particle limit by which 
we can introduce into general relativity a notion of a point-like particle with 
strong internal gravity without using Dirac delta distribution. Besides this 
limit, to deal with strong internal gravity we express the equation of motion 
^ ' in surface integral forms and calculate these integrals explicitly. As a result 

we obtain the equation of motion for a binary of compact bodies accurate 
through the second and half post-Newtonian (2.5 PN) order. This equation is 
derived in the harmonic coordinate. Our resulting equation perfectly agrees 
with Damour and Deruelle 2.5 PN equation of motion. Hence it is found that 
the 2.5 PN equation of motion is applicable to a relativistic compact binary. 
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PACS Number(s): 04.25.Nx, 



I. INTRODUCTION 

A relativistic compact binary system in an inspiralling phase is one of the most promising 
sources of gravitational waves for detectors under construction such as LIGO jjj], VIRGO 
@, GEO600 ||, and TAMA300 §. Among those detectors, TAMA300 has made rapid 
progress and has recently started observation ||. 

To detect gravitational wave signals directly, which are expected to be buried in noisy 
data streams, the matched filtering technique will be (and, in fact, has been) used. To use 
this technique, we have to construct reliable wave form templates, i.e., theoretical prediction 
of gravitational wave form accurate enough to match a data stream with templates without 
missing one phase of the expected signal. The equation of motion for the source system is 
one of the necessary ingredients for making such templates. 

In an inspiralling phase, since the separation of the binary is sufficiently larger than 
the radius of the stars, the interbody gravity is weak and orbital velocities of stars are 
small compared to the velocity of light. Hence the post-Newtonian approximation with a 
monopole-truncated object || description is useful. In fact, various authors have used it 
to derive an equation of motion of the system [^ -|27|,|7|-[14| . To obtain sufficiently accurate 



templates, we have to derive a highly accurate, say, the 4 PN equation of motion |]1^-|19| 
(though, the 2 PN templates [p0| |22f| might be enough to construct search templates [f2~3|1 ). 
What we are aiming at is to derive such an accurate post-Newtonian equation of motion. 
In this paper we derive the 2.5 PN equation of motion as a step toward this goal. 

The 2.5 PN equation of motion for two monopole-truncated bodies has been derived by 
various authors; Damour and Deruelle p4|-p7|1 used the post-Minkowskian formalism |28| . 



Kopejikin and Grishchuk worked on a spherical body J7[§|, Schafer ]ID| used the Hamil 



tonian approach, and Blanchet, Faye and Ponsot used the post- Newtonian approach IJI2" 
Their resulting equations of motion agree with each other. Then what can we add? It 
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seems to us that these derivations have some uncomfortable features. In the derivation in 
2^ , 2^ - ^7 , ID[ , I2| , though armed with the "Dominant Schwartshild" condition [Z7|, one needs 



to regularize divergences caused by their use of Dirac delta distribution. In the derivation 
free from any divergences one applies the post-Newtonian approximation even to the 
inside of a star (recall that a neutron star has strong internal gravity). Though the deriva- 
tion in [12| is mathematically systematic and well-devised, one assumes that a point mass 
follows a regularized geodesic equation. Instead, in our method we can take into account the 
strong gravity inside stars by use of the strong-field point particle limit [^9],[30l and surface 
integrals, without assuming a geodesic equation. 

One of our motivations to derive an higher post-Newtonian equation of motion with our 
method comes from the recent result on the 3 PN equation of motion, which has been derived 



by two groups in independent manners (see also [31]). Jaranowski and Schafer have used 



Hamiltonian approach and succeeded to determine the 3 PN order Hamiltonian [J 1] except 



two terms whose coefficients are arbitrary; the kinetic part ]TTj| and the static part [|33|]. The 
former has been recently determined by imposing poincare invariance on their equation of 



motion [3~3|. Blanchet and Faye have also succeeded to derive the 3 PN equation of motion 
IB|,[rjj using their generalized Hadamard partie finie regularization P5|P5 . Though their 



regularization method is declared to be well-defined unlike the method used by Jaranowski 
and Schafer JTT| , they found that there remains one arbitrary constant corresponding to the 
static part in the work of Jaranowski and Schafer. Although the method of Blanchet and 
Faye seems mathematically rigorous and well-devised, what is uncomfortable, besides the 
above arbitrary constant, comes from the fact that they rely on some "naturality" : They 
assume that a mass point follows a regularized geodesic equation, but the assumption should 
be verified. Hence, it is interesting to derive the 3 PN equation of motion using yet another 
method. Besides this, It is important to derive the equation of motion with various methods 
to confirm the equation, which appears through long calculations. 

In our previous paper ( |37| , referred as paper I henceforth), we have derived the 1 
PN equation of motion using the post-Newtonian approximation with the strong filed point 



particle limit p9| , |30| and the general equation of motion expressed in surface integrals |37| as 



in the Einstein- Infeld-Hoffmann approach |32| . The spin-orbit and quadrupole-orbit coupling 
forces which appear at the 2 PN order in our ordering have been also derived. Following 



the initial value formulation on the field by Schutz and Futamase |38|-f4T||, in this paper we 
shall derive the 2.5 PN equation of motion. 

This paper is organized as follows. In Sec. II, we will review the post-Newtonian approx- 
imation and the strong field point particle limit. In Sec. Ill we shall show how to solve the 
Einstein equation. There we shall define multipole moments of stars, and introduce a gen- 
eral form of the equation of motion. It is entirely expressed in surface integral form. From 
Sec. IV to Sec. VII, we shall derive the equation of motion for two compact bodies accurate 
though the 2.5 PN order, the order at which the radiation reaction effect first appears. It 
will be found that the temporal component of the equation of motion and its (functional) 
solution play an important role in our derivation. The section VIII is devoted to conclusion 
and discussion. Useful formulas and explanations on subtleties of our method will be given 
in appendices. 

Throughout we use units where c = 1 = G. x denotes Euclidean three vector. We raise 
or lower its index with Kronecker delta. 

II. POST-NEWTONIAN APPROXIMATION WITH STRONG FIELD POINT 

PARTICLE LIMIT 

In this section we shall review the strong field point particle limit associated with the 
post-Newtonian approximation. Also scalings of initial data for matter and a body zone will 
be explained. These scalings enable us to incorporate a self-gravitating point-like particle 
into post-Newtonian approximation and ensure that tidal effect is the lowest order effect 



which causes an internal motion of the stars except their spinning motion. See p9| , |30| for 
detail. 

In the inspiralling phase, the post-Newtonian approximation with a point particle de- 
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scription is suitable to derive the equation of motion. In the post-Newtonian approximation 
we assume the balance between Newtonian gravitational force and centrifugal force, i.e., 
m/L 2 ~ v 2 rb /L, where m, L, and 

Vorb £^1*6 typical scales of a mass of a star, the separation 
of the binary, and an orbital velocity of the star. Then a natural time scale r, which we 
call Newtonian dynamical time, is r = et |39(| , where e is a non-dimensional parameter and 
represents the smallness of the orbital velocity; 

■ doc doc , . 

^ = ~dt = e *F> (2 - 1} 

and we set dx l /dr = 0(e°). Therefore we take e as the post-Newtonian expansion parameter. 
Then, let e go to zero with r remaining constant. The meaning of r is that two events at 



the same (r, x l ) but having different e are in the nearly same phase in the orbit [39]. Then 
from the Newtonian force balance, m ~ e 2 when we keep the orbital separation constant. 

Now we achieve point particle limit by letting the radius of the star, say, R, shrink at 
the same rate as the mass of the star; m/R = 0(e°) = constant while R = 0(e 2 ) —>■ and 
m = 0(e 2 ) — > as e — ► 0. Because of finite m/R, we say this point particle has strong 
internal gravity. This scaling implies that the density changes proportional to e -4 (in the 
(£,£*) coordinate). The scalings of R and m motivate us to define the body zone of the 
star A (A=l,2), Ba = {x % \\x — za{t)\ < ^Ra} and a body zone coordinate of the star A, 
or A = e~ 2 {x % — z a (t)). z a {t) is a representative point of the star A, e.g., the center of the 
mass of the star A. R A is an arbitrary length scale (much smaller than L and is not the radius 
of the star) and constant (, i.e., dRA/dr = 0). With the body zone coordinate, the star does 
not shrink, while the boundary of the body zone goes to infinity. Then it is appropriate 
to define star's characteristic quantities such as a mass, a spin, and so on with the body 
zone coordinate. On the other hand the body zone serves us with a surface dB A , through 
which gravitational energy momentum flux flows and in turn it amounts to gravitational 
force exerting on the star A. Because the body zone boundary 8Ba is far away from the 
surface of the star A, we can evaluate the gravitational energy momentum flux at 8Ba with 
the post-Newtonian gravitational field. In fact we shall express our equation of motion in 



terms of integrals over OBa and be able to evaluate them explicitly. 

Let us move on to the scaling of the matter variables. As in the paper I, we consider an 
initial value problem to solve the Einstein equation. We take two nearly stationary solutions 
of the exact Einstein equation as the initial data for the matter variables (and gravitational 
field) . Then we assume that these solutions have the following scalings; the density changes 
proportional to e -4 , the velocity of spinning motion is 0(e). The slow spinning motion 
assumption is not crucial: In fact, it is straightforward to incorporate a rapidly spinning 
compact body into our formalism. The scaling of the density suggests that the natural 
dynamical time (free fall time) 77 inside the star may be rj = e~ 2 t. Then if we can not 
assume the nearly stationary condition on the stars, it is difficult to use the post-Newtonian 
approximation p9| , |30[ . 



From these initial data we have the following scalings of the star A's stress energy tensor 
components in the body zone coordinate, T^ T = 0(e~ 2 ), Tjf = 0(e~ 4 ), Tjf = 0(e~ 8 ). 
Here the underlined indices mean that for any tensor A 1 , A 1 = e~ 2 A l . In the paper I, we 
have transformed Tj^ u , the components of the stress energy tensor of the matter in the near 
zone coordinate, to using the transformation of Fermi normal coordinate at the 1 PN 
order. It is difficult, however, to construct the Fermi normal coordinate at an higher post- 
Newtonian order. Therefore we shall not use it. We simply assume that for Tj^ u (or rather 
A%, see Eq. Q), 

= 0(e" 2 ), (2.2) 
T^ = 0(e" 4 ), (2.3) 
T% = 0(e-% (2.4) 

as their leading scalings. We have found that these scalings are enough to derive the 2.5 PN 
equation of motion for two compact bodies. 

Henceforth we call the coordinate (r, x l ) the near zone coordinate. 
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III. FIELD EQUATION AND THE GENERAL FORM OF THE EQUATION OF 

MOTION 



A. Field equation 

As discussed in the previous section, we shall express our equation of motion in terms of 
surface integrals over the body zone boundary where it is assumed that the metric slightly 
deviates from the flat metric rj^ = diag(— e 2 , 1, 1, 1) (in the near zone coordinate). Thus we 
define a deviation field W as 

hr = rr - V=gg» u , (3.1) 

where g is the determinant of the metric. Our W" differs from the corresponding field in 
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m a sign. 

Now we choose the harmonic coordinate condition on the metric 

W\ v = 0, (3.2) 

where the comma denotes the partial derivative. Then, we recast the Einstein equation into 
the relaxed form, 

Uh^ = -16ttA^, (3.3) 

where □ = rj^d^d^ is the flat d'Alembertian and 

A" u = 6^ + x^ al3 ',*/?, (3-4) 
6^ = (-^)(T^ + tf L ), (3.5) 

16n 

Here, T^ v and t^ L denote the stress-energy tensor of the stars and the Landau-Lifshitz 
pseudo-tensor E2[. yV-vap originates from our use of the flat d'Alembertian instead of the 



curved space d'Alembertian. In consistency with the harmonic condition, the conservation 
law is expressed as 



A>% = 0. (3.7) 



Now we rewrite the relaxed Einstein equation into an integral form, 

3 A^(r -e\x-y\,y k ;e) 



h» u = 4 d 6 y v ,J *"* ' (3.8) 

JC(T,x k ;e) \X — y\ 

where C(r, x k ;e) means the past light cone emanating from the event (r, x k ) (about the 
dependence of C on e, see ||3~9| ). We ignored the homogeneous solution, which corresponds 
to the initial data for the field, for simplicity. Even if we take random initial data for the 
field |58| supposed to be of the 1 PN order |40| , they are irrelevant to the dynamics of the 
binary system up to the radiation reaction order inclusively [ff(] . 

We shall solve the relaxed Einstein equation as follows. First, we make retarded expan- 
sion of Eq. ( |3.8|) and change the integral region C to iV (, N denoting near zone), a r = 
constant spatial hypersurface (note that such retarded expansion is understood to be formal, 
since it gives divergent integrals. Because such divergent integrals do not appear up to the 
2.5 PN order and we are concerned with the 2.5 PN equation of motion in this paper, we 
use formal retarded expansion. See below). 

= 4E ^"j^Alf-^A^y^e). (3.9) 

Here we added a subscript N to A M ^ to clarify that here are components in the near 
zone coordinate. Simple retarded expansion gives divergent integrals at the 3 PN order. In 
this paper we take the near zone as a sphere centered at some fixed point and enclosing the 
binary system. The radius of this sphere is set to be IZ/e, where TZ is arbitrary but much 
larger than the size of the binary. Then we find that up to the 2.5 PN order, the terms 
dependent on TZ take forms such as 7Z~ n , d T \nlZ, <9jln7£, d T TZ n , or dilZ n (n > 1), where 
8^ = d/dx^. Hence we simply let 1Z go to infinity. 

Second we split the integration into two parts; contribution from the body zone Ba, and 
from elsewhere, N/B. Schematically we evaluate the following type integration (we omit 
indices of the field), 
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h = h B + h N/B , (3.10) 



^ = e 2^ / rf a A — 2 ^ M_ ra , 3.11 



h N/B =[ A ,/ (r f n , (3.12) 
Jn/b \x — y] 1 n 

where r A = x — ^4. We shall deal with these two contributions successively. 

1. Body zone contribution 

As for the body zone contribution, we make multipole expansion with care over the 
scaling of the integrand, i.e., in the body zone. For example, n = part in Eq. ( |3.9|) , 

(pr r>k k ofkl k I \ 

^ + e 2 ^ + e 4 ^f^M + 0(e 10 ), (3.13) 
r A r A 2r A J 

/ pi jki r k \ 

hs n = =^ £ ^ + e 2 ^ +0(e 8 ), (3.14) 



A=l,2 



h%n,o = 4, 2 E f 1 + + + <*"). (3.15) 

A=l,2 V r ^ r A Ar A J 

Here the quantity with a hat is symmetric and tracefree on its dummy indices, and r A = \r A \. 
To derive the 2.5 PN equation of motion, h TT up to 0(e 9 ) and /i T! and h % i up to 0(e 7 ) are 
needed. 

In the above equations we defined multipole moments as 

l L A = e 2 / <Pa A XZa%, (3.16) 
Jb a 

J% = e 4 / d 3 a A A^c4 (3.17) 

7^ii 8 / j3„, A 4? ^ 



Z^=e 8 d 3 a A A%a% (3.18) 

where we introduced multi index notations L = iii 2 • • • H and a A = cc A cc A ■ ■ ■ or A . Then 
P A = I A , D\ = P A , P A = J\. We simply call P A the four momentum of the star A, P\ the 
momentum, and P A the energy. Also we call D k A the dipole moment of the star A, and I A 
the quadrupole moment of the star A. 



Then we transform these moments into more convenient forms. By the conservation law 
Eq. ( |3.7| ), we have 



A 



N 



a| = (aM), t + (aM), + ^ 



N ) 



(3.19) 
(3.20) 



where v\ = z l A , an overdot denotes r time derivative, and y A = y — z A . Using these equations 
and noticing that the body zone remains unchanged (in the near zone coordinate), i.e., 



0, we have 



7«J 



P l A = P T A v\ + Q\ + e 



2 dD\ 
dr 



J A 



+ v\D 3 l + -e 



2 Q% 



;dD 



1 d 2 T lJ 

eP A v A v A + 2 e dr2 + 2ev A ^ 



j) 



.1 A^ V A n i) 



l TV 

dr UA 



+ e<Q^ + e<R^> + -e ^ 



7 kij 
^A 



-/1 A — /i^ , 



(3.21) 
(3.22) 



(3.23) 
(3.24) 



where 



A/J = 2e 4 / rf 3 a A a^A- w . 



JdB A v 7 



= e" 4 / (Aft - «* A]J) j^ri, 



(3.25) 
(3.26) 
(3.27) 



and 



A A - e J A V A + e V A J A + K A + e —, • 



(3.28) 



| ] and ( ) attached to indices denote anti-symmetrization and symmetrization. M A is the 
spin of the star A and Eq. (|3.21|) is the momentum-velocity relation. In general, we have 



jLi 
J A 



L A 



( W ) 21 T (L-l[i t )i] 



l + l 



r{Li)j 



+ 



21 

7+T 



r(L-l[ii)i]j 



+ z 



(Lj)i 



+ 



21 

7+T 



,(L-l[ii)j}i 



(3.29) 
(3.30) 
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and 



z (Li)j + z (Lj)i = e 2 v pL)j + e 2 y ' + _±_ e 4^A_ + e -^+2 jR M^) ) (3.32) 

where / is the number of indices in the multi index L. 

Now, from the above equations, especially Eq. ( |3.23| ), we find that the body zone con- 
tribution, h^ n=Q) are of order 0(e 4 ). Note that if we can not or do not assume the (nearly) 



stationarity of the initial data for the stars, then, instead of Eq. ( |3.23| ) we have 

where we used the dynamical time 77 (see Sec |T|). In this case the lowest metric differs from 
the Newtonian form. From our (nearly) stationary assumption the remaining motion inside 
a star, except the spinning motion, is caused only by the tidal effect by the companion star 
and from Eq. ( [3.231 ), it appears at the 3 PN order [[30 . 



To obtain the lowest order h^ n=Q , we have to evaluate the surface integrals Q^,R^. 
Generally, in h^ n=Q the moments J^ z and appear formally at the order e 2Z + 4 and e 2l+2 . 
Thus Qa and R L ^ appear as 

r L/)Li 

n Bn=o ~ r e 2l+1 H , 

r A 

n Bn= $ ~ he 2^+i H ' 

where we omitted irrelevant terms and numerical coefficients. Thus one may expect that 
Qa and appear at the order e 4 for any L and we have to calculate an infinite number of 
moments. This is not the case. We address this problem in the appendix [X]. The important 
thing here is that e 4 Q^ and e^R^ are at most 0(e 4 ) in h^ n=Q . 

Finally, since the order of h^ n {n > 1) is higher than that of h^ n=0 , hf£ = 0(e 4 ). 
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2. N/B contribution 

About the N/B contribution, since the integrand A^T = —gt^ L + x^ ualS ,a/3 is at most 
quadratic in the small deviation field W, we make the post-Newtonian expansion in the 
integrand. Then, basically, with the help of a potential g(x) which satisfies Ag(x) = f(x), 
A denoting Laplacian, we have for each integral (, e.g., n = term in Eq. ( |3 . 1 3 ) ) 



/ d 3 y-l^- = -4ng(x) + / dS k 

J N/B \X — y\ Jd(N/B) 



1 dg(y) d ( 1 



(3.33) 



\x — y\ dy k dy k \\x — y\ t 

The last surface integral really contributes to the 2.5 PN gravitational field. For n > 1 terms 
in Eq. ([3.12 ), we use potentials many times to convert all the volume integrals into surface 



integrals and "— 47rg(a ? )" terms. 

Now the lowest order integrands can be evaluated with the body zone contribution hp/, 
and since is 0(e 4 ), we find (see appendix ||) 

An = 0(e 6 ), (3.34) 
A" = 0(e 6 ), (3.35) 
K = (~9)4l + 0(e 8 ) = e 4 ^- (5\5h - l^fo) + 0(e 5 ). (3.36) 

where we expanded h% in an e series; 

n=0 

Similarly, in the following we expand h^ u in an e series. From these equations we find that 
the deviation field in N/B, h^, is 0(e 4 ) (It should be noted that in the body zone is 
assumed to be of order unity and within our method we can not calculate there explicitly. 
To obtain in the body zone, we have to solve the internal problem). 

The Landau-Lifshitz pseudo- Tensor expanded with h^ v is listed in the appendix M. 



B. General Form of the Equation of Motion 

From the definition of the four momentum 
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PZ(r) = e 2 [ d 3 a A A T J, (3.37) 

JB A 

and the conservation law, Eq. (|3.7|) , we have the evolution equation for the four momentum; 

^A = - e -±t dS k K^ + t- A v k A i dS k A T £. (3.38) 
dr JdB A JdB A 

Here we used the fact that the size and the shape of the body zone are defined to be fixed 
(in the near zone coordinate). 

Inserting the momentum- velocity relation, Eq. ( |3.21| ), into the spatial component of Eq. 
( |3.38| ), we obtain the general form of the equation of motion for the star A; 



dr JdB A JdB A 

dSkA 1 ^ 

dB A JdB A 



dQ\ ^fD^ (33g) 



dr dr 2 



All the right hand side terms in Eq. ( |3.39p except the dipole moment are surface integrals. 



We can specify the value of D A freely to determine the representative point z a (t) of the 
star A. In this paper we take D\ = and simply call z\ the center of the mass of the star 
A. Note that to obtain the spin-orbit coupling force in the same form as the previous works 
f^j43"lj44j , we have to make another choice for z\ (paper I). For completeness we show in 
the appendix |C] the spin-orbit coupling force on the condition that D l A = 0. Henceforth we 
shall restrict our attention on the mass monopoles. 

In Eq. ( |3.39| ), P\ rather than the mass of the star A appears. Hence we have to derive a 
relation between the mass and P r A . We shall achieve this by solving the temporal component 
of the evolution equation ( |3.38| ) functionally. 

Then since all the equations are expressed with surface integrals except D l A to be spec- 
ified, we can derive the equation of motion for the strongly self-gravitating star using the 
post-Newtonian approximation at least before the order where we have to be concerned with 
the internal problem. 
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C. On the Arbitrary Constant Ra 



Since we introduce the body zone by hand, the arbitrary defined body zone size Ra 
seems to appear in the metric, the multipole moments of the stars, and the equation of 
motion. More specifically Ra appears in them because of 1) the splitting of the deviation 
field into two parts (, i.e., B and N/B contributions), the definition of the moments, and 2) 
the surface integrals that we evaluate to derive the equation of motion. 



B and N/B contributions to the field depend on the body zone boundary, €Ra- But W 
itself does not depend on cRa- Thus it is natural to expect that there are renormalized mul- 
tipole moments which are independent of Ra since we use non-singular matter sources. One 
possible practical obstacle for this expectation might be log(ei?^) dependence of multipole 
moments. We have checked that at least up to the 2.5 PN order, there is no such log-term. 

Though we use the same symbol for the moments henceforth as before for notational 
simplicity, it should be understood that they are the renormalized ones (we use the symbol 
"Pa" for the renormalized P4). 



Since we make integrations over the body zone boundary, in general the resulting equation 
of motion seems to depend on the size of the body zone boundary, €Ra- But it does not 
depend on Ra- 



1. Ra Dependence Of The Field 



2. Ra Dependence Of The Equation Of Motion 



In the derivation of Eq. (|3.39|) , if we did not use the conservation law (|3.7|) until the 



final step, we have 
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dQ t A+f ?d'D\ 



dr dr 2 

= e" 4 J d 3 y^, u - e~\\ jf d 3 yA T N \ u + e" 4 ^ #vA%*vty • (3-40) 

Now the conservation law is satisfied for whatever value we take as Ra, then the right hand 
side of the above equation is zero for any Ra- Hence the equation of motion Eq. (|3.39|) does 
not depend on R A (Similar argument can be found in [p2| ). 

Along the same line, the momentum-velocity relation ( |3.21| ) does not depend on R A . 

Up to the 1 PN order in paper I, we have explicitly shown the irrelevance of the equation 
of motion to Ra by checking the cancellation among the Ra dependent terms. 

IV. NEWTONIAN AND THE 1 PN EQUATION OF MOTION 

In this section we shall list the Newtonian and the 1 PN equation of motion for later 
convenience. As for the details of the derivation in our method, see the paper I. 
Let us review how to derive the field and the equation of motion order by order. 
First of all, from Eqs. (|3.34j) , (|3.35| ), and the time component of Eq. ( |3.38|) , we have 



iS.o^. (4.1) 

Then we define the mass of the star A as the integrating constant of this equation; 

m A = \imP T A . (4.2) 

m A is the ADM mass that the star A had if A were isolated (we took e zero limit in Eq. 
Q4.2D to ensure that the mass defined above does not include the effect of the companion 
star and the orbital motion of the star itself. By this limit we ensure that this mass is the 
integrating constant of Eq. (|4.1|) . Some subtleties about this definition are discussed in the 
appendix 0). By definition ttla is constant. Then we obtain the lowest order h TT \ 

ih r = aJ2—- (4-3) 

A=l,2 TA 
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Second, from Eq. (|3.26| ) with Eqs. ( |3.34|) and ( |3.35|) we obtain Q l A = at the lowest 
order. Thus we have the Newtonian momentum- velocity relation P A = rriAV A + 0(e 2 ) from 
Eq. ( |3.21| ) (we set D\ = 0). Then from Eq. ( |3.39| ), we can derive the Newtonian equation 
of motion. This completes the Newtonian order calculations. 

Next, at the 1 PN order, we need eh TT , 4/1^". The n = 1 term in the retardation expansion 
series of h TT , Eq. ( |3.9|) , gives no contribution at the 1 PN order by the constancy of the 
mass TiiA, i.e., ^h TT = 0. 

Now we obtain Ji Tl from the Newtonian momentum-velocity relation. 

A=l,2 T A 

From 4h TU and Eqs. (|B3| ) and ( |B4l) we can evaluate surface integrals in the evolution 
equation for P r A at the 1 PN order. The result is (for the star 1) 

dP[ -^4(U+—)+0(e% (4.4) 



dr dr \2 r 12 

where fu = z\ — z 2 and r\ 2 = \r\2\, and we used the Newtonian equation of motion. From 
this equation, we have the mass-energy relation at the 1 PN order, 



PI = mi 



, 2 / 1 2 , 3m 2 
1 + e ~v{ 1 



0(e 3 ). (4.5) 



,2 1 r 12 

Then we have to calculate the 1 PN order Q A from Eqs. (p3|) , (|B4|) and (|3.26 ). The 
result is that Q\ = e 2 m A v l A / (6€Ra)- As Q\ depends on R A , we ignore it (see Sec. |III Q) . As 
a result we obtain the momentum- velocity relation at the 1 PN order, P\ = P A v A + 0(e 3 ) 
from Eq. ( EHHD . 

Now as for 4 /i y , we first calculate the surface integrals Q l ^R l ^R k A ^ from Eqs. (|3.26| ) 
and ( |3.27| ) using Eqs. (|H^D, (P^|), and ( |E4] ). We then find that they depend on R A , hence 
we ignore them and obtain 

h% = Ae A £ ^>i + (e 5 ). 

A=l,2 T 'A 

To obtain s h TT and 4/1^, we have to evaluate non-compact support integrals for 6^]v/b 
and 4h l ^ B , and n = 2 term in Eq. (|3.9| ) for h TT . It can be done with the help of Laplacian 
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inverse; Alnri = 1/rf and AlnS* = \jr\r 2 [p5 ]. Using the Newtonian equation of motion 
we obtain 



If 



pr 

4e 4 £ £A 

A=l,2 r ^ 



. , 171 A r/^ ^ \2 oi 17111712 _ _ 

-2 > {(ru-tu) -u A } + 2 — 2 — n 12 -{n 1 -n 2 ) 

A=l,2 r ^ r 12 



+ 7E 



+ 14 



17li17l2 



A=l,2 T A 



r x r 2 



- 14 



17111712 



E- 

r 12 AtT,2^ 



(4.6) 



fcij = 4e 4 £ 



m A v A v A 



A=l,2 T A 
2 



v - 4 j »mim 2 

mi7T7.2 



A=l,2 r A 



2 J 

n 12 ra 12 



g 2 



n l 2-n l ) {k (n l 2 + n 2 ) l) 



(4.7) 



Here 5 = r x + r 2 + r 12 , n 12 = f 12 /r 12 , n A = f A /r A . 

Inserting the field 4W, §h TT into Eqs. (|B4j ) and ( |B5|) and evaluating the surface integrals 
in Eq. ( |3.39| ), we obtain the 1 PN equation of motion; 



TO! 



dv\ 
d7 



17li17l 2 i 
2 n 12 

n 2 

2 m 1 m 2 
r\ 2 



9 ^9 3,^ _ w _ 5mi 4m 2 

-v? - 2i; 2 2 + -(n 12 ■ v 2 ) 2 + 4(^ ■ ?T 2 + + 

2 r 12 r X2 



+ r (4(n 12 • v x ) - 3(n 12 • v 2 )) 



(4.8) 



where we defined the relative velocity as V = v\ — v 2 and we used the Newtonian equation 
of motion and Eq. (^.5|) . 

Finally let us give a summary of our procedure. With the n — 1 PN order equation of 
motion and h) iV in hand, we first derive the n PN evolution equation for P A . Then we solve 
it functionally and obtain the mass-energy relation at the n PN order. Next we calculate Q\ 
at the n PN order and derive the momentum-velocity relation at the n PN order. Then we 
calculate Q% and R L A ij . With the n PN mass-energy relation, the n PN momentum-velocity 
relation, Q A , and R A ^, we next derive the n PN deviation field h^ v . Finally we evaluate 
the surface integrals which appear in the right hand side of Eq. ( |3.39|) and obtain the n PN 
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equation of motion. In the above calculations we use n-1 PN order equation of motion to 
reduce the order, e.g., when we meet e 2 dv\/dr in the right hand side of the equation motion 
and we have to evaluate this up to e 2 , then we use the Newtonian equation of motion and 
replace it by — e 2 m2r l l2 /rf 2 - Basically we shall derive the 2.5 PN equation of motion with 
the procedure as described above. 



V. THE 1.5 PN EQUATION OF MOTION 

In the appendix we show that Q L A % and do not contribute to W up to the 2.5 
PN order. There, we also derive the momentum- velocity relation up to the 1.5 PN order 



P l A = P T A v\ + 0(e 4 ), 



and the relations between and some moments; the results are that Z l j[ = e 2 P\v % A v A + 
0(e 4 ) and Z A j = e 2 Mfv j } + 0(e 4 ). Henceforth we shall use Eqs. O and ( |AT"0l ). 



Now let us derive the 1.5 PN equation of motion, which is 0(e 3 ) correction to the 
Newtonian equation of motion. Thus as the integrands in Eq. (|3.38|) , we shall use r[— fftLi] 
and hence we have to obtain jh TT and ^h^ 1 . The fields that we shall evaluate are 

d r ,„ t 2 „ d 3 



(5.1) 

/ l « = 4< / l "_4e 5 i dS M {-gq L ]+0{^), (5.2) 

JdN 

h ij = 4<^ r - 4e 5 -f- f £ P t aV>'a + e 4 f d^[-gt[ L ]\ + 0(e% (5.3) 
dr ^ x ti, 2 Jn / b J 

Here < n h^ u denotes the field complete up to the order e n . 

Now let us follow the procedure described in the previous section. From Eqs. ( [B6D and 

(B7|), the mass-energy relation at the 1.5 PN order becomes 

P T A = (P T A )<iPN + 0(e i ) 

Here {P£)<ipn stands for the 1 PN order mass-energy relation, Eq. ( |4.5| ) (without 0(e 3 ) 
symbol). 



A straightforward calculation for the N/ B contribution results in 

Pa 



^ E 



A=\,2 



TA 



d 2 



2 E Q-5 tfM 

A=l,2 UT 



T3T T3T 
7 V" a b 

0,6=1,2 raTb 
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prpr I 
r 12 ^,2^ 



<5 



3 dr 
4e 4 



V A=1,2 

Plv\ 



A=l,2 Ta 



< 5 ^' = 4e 4 £ 

.4=1,2 



+ e 4 



pr2 

E 1 "' " ' ! 



8P[P^ 



A=l,2 'A 



1 J 

2 n A n A „ ~ c~ n 12 ra 12 



_ 9.5(3) jij 
ze J or6- 



P{Pl 



S 2 



(n 1 2-n i y k (n 1 2 + n 2 ) l) 



(5.4) 
(5.5) 



(5.6) 



Here I% b = £a=i,2^144 and (")/ denotes d n I/dT n . 

Note that 5/t u and 7/i TT depend only on time, thus from Eq. 
at the 1.5 PN order becomes 



dv\ 
~dT 



dv\ 
dr 



+ 0(6 4 ). 



the equation of motion 



<1PN 



VI. THE 2 PN EQUATION OF MOTION 



In this and next section, we shall derive the 2.5 PN equation of motion using the proce- 
dure described in the section [IV]. The main problem that we have to solve is to derive h^, B . 
Thus, from Eq. ( |3.33|) , the problem is reduced to solve the Poisson equations Ag(x) = f(x). 
Fortunately, all the solutions of the Poisson equations that we need to derive the 2.5 PN 



equation of motion are obtained in [12.11 
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As the first step of the 2 PN order calculation, we have to derive §h T% which will be 
needed to obtain the 2 PN order evolution equation for P T A . 



<«h Ti = 4e 4 E 



Plv 



A v A 



+ 4e b 



d 3 y 



e[-gt 



Tl 1 

llJ 



2e fc 



Q2 



E pT AV l A r A 



V A=1,2 



A=i,2 r ^ J n/b ' \x-y\ 

Now we shall show calculations for the contribution in detail. The integrand is 



(6.1) 



■lGngt2 L ] = 2;// 7 '-' - 'pi 



Then the first integrand contributes 
2 



167T J N/B 



3 4^-' ,fc4 

d y- 



x-y\ 



h r[k,i] 2 

7T 



A=l,2 



pr2 [fc 
-'A U A 



- APlP^d 2A d lk \nS - APlP;vfd^d 2k \nS 



(6.2) 



where d A i = djdz\. In the second equality of the above equation, we used the following 
formulas, 



d 3 y r* 



A' A 



n/b \x — y\ r\ 



n/b \x — y\ 



(6.3) 



/ 

JN 



d 3 y r\r 2 



d 3 y- 



d d 



(AlnS), 



I n/b \x — y \ r\r 2 J n/b * \x — y\ dz\ dz 2 
and Eq. (|3.33|) . Applying the same formulas to the other integrand, we find at last, 

" -gtli] ^ Pa 2 r,* .-t , *..n , 4p i rp 2 r , 



(6.4) 



4 / ; dY- 

Jn/b \x — y\ 



E ~± {(^ ■ v A )n\ + 7v\} + -^{vik + v 2k ) (5 ki - n\ 2 n k 2 ) 



A=l,2 'A 



+8P[P 2 r (v[ +vi)[-— E — + — 



•16 



+ 12 



pin 

s 2 

f_l f_2 

S 2 



{vf(n* 2 - n\)(n k 2 + n k 2 ) + t; 2 fc « - n k ){n\ 2 + n 2 )} 
{^K 2 - n k )(n\ 2 + n 2 ) + v k (n\ 2 - n\)(n k 2 + n k )) (6.5) 



We show e h Tl in the appendix ||. 
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A. The 2 PN evolution equation for P A 

In this subsection we shall derive the 2 PN evolution equation for P A . The surface 
integrals we have to evaluate are 

A pr 



_ ; / dS k6 [-gt T L k L ] + e 2 v k A I dS k6 [-gt T LLl 
/ dS k8 A k N T + e\ k A £ dS k8 A T N T + 0(e 5 ). (6.6) 



e 4 



/SBa J dB A 

Here 0(e 2 ) terms become 

-/ d^ 6 [-^li]+«f/ ^[-^ll] = [4(nia • #i) - 3(n 12 • tf 2 )] . (6.7) 

JdBi JdBi T^2 

This is the 1 PN order evolution equation with the mass tua replaced by the energy P A . 

From the 2 PN order, XaT comes into play. We can see this by the definition of P A , 
Eq. ( |3.37| ). Noticing Xn ,ap — Xn %3 we rewrite Eq. ( |3.37| ) as 



PJ(r)=e 2 / d 3 a A Q}J + e- 4 i dS iX T \r (6-8) 

JB A JdB A 

Since Xn ,<*P * s °f or der e 8 , the surface integral of it first appears at the order e 4 , i.e., at 
the 2 PN order. 

Now since 9^ and Xjv"^ \ap are conserved separately, i.e., O^f tfl = = Xn"*^ ,aPn, we 
shall deal with 0-^ and Xjv separately. Thus we split P A into two parts; P Ae and P\ x - 

PUr)=e 2 f rfVeJ, (6.9) 

J Ba 

P^(r) = e 2 f d 3 a AX T P , a p- (6-10) 

J B A 

We consider the evolution equation of P A in the appendix |]. The evolution equation for 
P Ae becomes 



<1.5P7V 



dPX e _ dPJ 
dr \ dr 

- 6 4 / dS k8 [-gt k L \] + e 4 v k A I dS k8 [-gtll]. (6.11) 

JdB A JdB A 

From Eqs. ( |B9"|) and (|B10| ), the relevant combinations of h^ v are; 
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4 h Ti > j 4 h kl > m , 4 h TT ^ 4 h jk , T , 6 h Ti > j 4 h TT > k , Q h TT '\h Tj ' k , 4 h TT 4 h TT ^ 4 h Tj ' k 
4 h Ti 4 h TT > j 4 h TT > k , 6 h TT ^ 4 h TT ' j , 4 h Ti ' j 4 h Tk ' 1 , 4 h TT ^ 4 h jk ' 1 , 4 h TT 4 h TT ^ 4 h TTd . 

Evaluating surface integrals of these ten types of integrands, we obtain, 

mim 2 



9Bi 



dS k8 [-gt 



and 



rk 1 
LLl 



9(^12 • V2) 3 5(ni2 • vi)(n 12 • v 2 f 4vf(n 12 ■ Vi) 



2 2 

4(^i • v 2 )(n 12 • vi) 
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+ 3(^i • v 2 )(n 12 ■ v 2 ) 

7vl(n 12 -v 2 ) 7vl(n 12 ■ vi) 



+ 



r\ 2 
m\m 2 
r\ 2 



-2(n 12 ■ v 2 ) - 



2(n 12 • vi) 



19(ni2 • v 2 ) _ 31(ni2 • gO _ 7(n u ■ V) 
2 6 2 



(6.12) 



9Bi 



mim 2 



r 



12 



-«i(ni2 • #2) + 



7(ni2 • tfi)(ni2 • W2) 2 . 4:vf(n 12 -vi) 



+ 



+ 



+ 



(vi ■ v 2 )(n 12 -v 2 ) + 

Amiml(ni 2 ■ Vi) 
3 

28mlm 2 (ni 2 ■ vi) 



2 15 

8(^1 • v 2 )(n 12 ■ vi) 5vi(n 12 -vi} 



(6.13) 



Combing these results, the evolution equation for P T A becomes 



'dR 



ig 



dr 



—e 



2 m\m 2 



<2PN 



+ e 



4 m 1 m 2 
r\ 2 



[4(n i2 • vi) - 3(n i2 • v 2 )\ 



2 2 
+ 4(tTi • v 2 ){n 12 ■ V) + 5vl(n 12 ■ v 2 ) - Avl(n 12 ■ vi) 



+ e 



4 mim2 



^12 



+ e 



4 m\m 2 
'12 



[-10(ni2 • ui) + 11(^12-^2)] 
[-4(n i2 • v 2 ) + 6(ni2 • #i)]. 



(6.14) 



We can functionally integrate this equation and the result is 
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P[ e = m l (l + e 2 2 T l + e\T 1 ) + 0(e 5 ). 



(6.15) 




The evaluation of the first three terms may call for great care. Let us show the results. 

- 32P 1 i P 2 j + <9 2 <<9 2 ;) In 5 

32(A • P 2 ) „ 1 . . 

+ 1 I 2J £ (6-22) 

dr 12 A=1.2 r ^ 



where the last term comes from the boundary term (surface integral term) in Eq. (|3.33|) . 

-^—,h l hh TT « = 4^ (>' - 3n\n{) 
n/b Att\x — y\ r{ v ' 

+ \QPlPlv\v{d 2 id^\nS 



4 pr3 

37? 



1QP[P;V! , ATD r2 V r( 4 1 



+ 4P^PJ -—j. + — — + (1 «- 2), (6.23) 
3 \ 3r!rf 2 6r 2 rf 2 / 

where (1 <-» 2) means that there are the same terms but with the labels 1 and 2 exchanged. 
The last two terms arise from the body zone boundary term in the Eq. ( |3.33 ). H x and K x 



have been derived by Blanchet, Damour, and Iyer ||46|| . They satisfy the following poisson 
equations in the sense of usual function in N/B; 

AK, = 2didj (1) dftQnn), (6.24) 
AHi = 2didj (—) dpd^ In S. (6.25) 



Their explicit expressions are 



1 1 1 r 2 rf 2 r\ 

12 

3i 

Hi = ~ «-> 3 ~~ 73 7l2l r>„2„2 + oX ^ 73l3~ + r. 3„2 o„3„3 ' (6.27) 



^1 — 3 H 2 2 1" o 2 2 + o 2 3 + o 3 2 ' (6.26) 

r 2 r 2 r i2 r i r 2 *r(r{ 2 2rfr 2 2r 2 r{ 2 



111 r 2 r 2 3r| 

2rf 4r? 2 4rfr i2 2r\r\ 2 2rfr i2 4rfr? 2 ' 2r\r\ 2 2r\r\ 2 



We note that they go to zero when r — > oo. The third term of Eq. ( |6.21| ) can be evaluated 

as 

-d 3 y 



n/b Att\x — y 



./' ;v ;/'",, 2 £ PI 2 ^>i(3^i^-^'A A )ln^ 



A=l,2 



+ 8 £ Pl 2 a^ fc In r A 

A=l,2 
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3r 12 A tt,2 r A 

of the star A, 



+ mPlP^vy^kd^ + v%v l 2 d 2k d 2l ) \nS 
+ 16P[P 2 T (a'ld 1 k + a k 2 d 2k ) In S 



a A — ~r~ 

A dr 



The last term again comes from the boundary term in Eq. (|3.33 ). 
The complete h^ u are shown in the appendix |]. 

C. The 2 PN equation of motion 

To derive the 2 PN equation of motion, we have to derive the momentum-velocity relation 
valid up to 0(e 4 ). Then, when the integrand is gA^, Q\ becomes 

Q\ = Qall + Qa x (6-29) 
• _ m\m 2 v\ _ m\v\v\ _ 2m\m 2 v % 2 
^ lLL ~ 3eR iri2 ~ 15eRi 3ei?ir 12 ' ^ ' ' 

Qa x = (6.31) 

where Qall an d Q\ x come fr° m the integrands —gt^ L and x*n ,<*p- Thus the momentum- 
velocity relation becomes 

ne = « + 0(e 5 ), (6.32) 

P\ x = P T Ax v\ + 0(e 5 ). (6.33) 

Now, by the divergencefree nature of Xjv ,a/3, Pa x is itself conserved (see the appendix 
|F]). Hence we shall consider P%q only. Up to the 2 PN order, we have to evaluate 
dP i 2 r 2 r 

-T = -E^ dS M+2n [-gtf L ] + J2 t 2n v A f dS M+2n [-gt* L ] + 0(e 5 ). (6.34) 

uT n=0 JdB A n=1 JdB A 

Evaluating the above surface integrals and using the mass-energy relation valid up to the 
2 PN order, we obtain the 2 PN equation of motion. We show it and the 2.5 PN order 
correction in the next section. 
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VII. THE 2.5 PN EQUATION OF MOTION 



In this section we shall derive the 2.5 PN equation of motion. 

Now we first determine the mass-energy relation up to the 2.5 PN order. The surface 
integrals we have to evaluate are 

dP Ae ( dP\e \ 



dr 



dr 



<2PN 



e 5 / dS k9 [-gt k L T L } + e 5 v k A <f dS k9 [-gtZ]. 



(7.1) 



ldB A JdB A 

The evolution equation for PT is shown in the appendix [F]. Evaluating these surface inte- 
grals, we find that 

0. 



AO 



dr 



2.5PN 



Thus we have the mass-energy relation (of P Ae part) 



Pie = mi l + t 2 2 T 1 + t\T 1 +0(e 6 ). 



(7.2) 



Here the definition of n Y A are given by Eqs. ( |6.16|) and (|6.17|) . We discuss an interesting 
interpretation of this mass-energy relation in the appendix ||. 

Next we have to derive the momentum-velocity relation. By evaluating explicitly, we 
find that Q A = 0. Hence we have 

P\ = P T A v\ + 0{e & ). 



Finally we have to derive the deviation field h^ u up to the 2.5 PN order; gh TT , •jh}" J . We 
have 



e 7 7 h TT + e 9 9 h TT 



4 7 d 
3 dr 



E nv\ - 



A=l,2 



pr pi 
f_l r 2 

r\2 



A=l,2 T A 



n l A n J A - -cf J 



E m ^ A 



JL 9 JUL 

+ (irrelevant terms), 



(7.3) 
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7 uri 



e' 7 h 



e\h\ + e 7 7 h\ 



—e 



,2&_ 
3^ 



(7.4) 



A=l,2 



-2e 5(3) / z 



6 35r3 A=l,2 

+ (irrelevant terms). 



x - 2 "^1^2 \ - 

L ™^a?A - — — ^ 



2r 12 



A=l,2 



(7.5) 



Evaluating the surface integrals in Eq. ( |3.39p with the deviation field above, we at last 
obtain the 2.5 PN order equation motion. 



mi 



dv\ 
dr 



m x m 2 i 
+e — 5— n\ 2 



12 



-v\ - 2v\ + 4(^i • v 2 ) + - (n 12 ■ v 2 f + + 
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+e 
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15 
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^i(ni 2 ■ ^2) + 4w 2 (^i2 ■ Vi) - 5v 2 (n 12 ■ v 2 ) - A(vi ■ v 2 ){n 12 ■ V) 

9 
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(7.6) 



This is the 2.5 PN equation of motion for the relativistic compact binary. We used the 
harmonic coordinate to derive this equation. This equation perfectly agrees with the previous 
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works [|4],[^,p,|l^] . Thus it is found that up to the 2.5 PN order the post-Newtonian equation 
of motion is applicable to the binary system even whose constituent stars have strong internal 
gravity. 

VIII. CONCLUSION AND DISCUSSION 

In this paper, we derived the equation of motion for two compact bodies accurate through 
the 2.5 post-Newtonian order, where the radiation reaction effect first appears. In the paper 
I, we have derived the spin-orbit and quadrupole-orbit coupling forces which are the only 
multipole-orbit coupling forces up to the 2.5 PN order in our ordering. Combing these 
multipole-orbit coupling forces and the result in this paper, we obtain an equation of motion 
appropriate to an inspiralling binary star system up to the 2.5 PN order. We have derived 
all these accelerations in the harmonic coordinate. 

To deal with a strongly self-gravitating object such as a neutron star, we used the gen- 
eral form of the equation of motion with the strong filed point particle limit. We imposed 
the scalings on matter and field variables on the initial hypersurface. This scalings enabled 
us to introduce a self-gravitating point-particle and ensured the (quasi-)stationarity of the 
stars. And in turn, the (quasi-)stationarity ensured that the internal motion of the star is 
stimulated primary by the tidal effect. In our method, we did not a priori assume a geodesic 
equation, but used only the local conservation law of the total stress energy-momentum ten- 
sor (matter's stress-energy tensor plus Landau-Lifshitz stress-energy pseudo-tensor). Using 
this conservation law, we introduced the general form of the equation of motion which is ex- 
pressed entirely in surface integrals. The evaluation of these integrals were made explicitly. 
The resulting 2.5 PN equation of motion agrees with Damour and Deruelle 2.5 PN equation 



of motion p^ , |25| j8|JT2]] . Hence we find that the 2.5 PN equation of motion is applicable to 
a relativistic compact binary, modulo the scalings imposed initially on the matter and field 
variables. 

Throughout our derivation, the evolution equation for the time component of the four 
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momentum of the star A, PJ, has played an important role. As a by-product, we obtain an 
interesting and natural relation between P T A and the binary's characteristics, Eq. (|H3| ) . 

We ignored the dependence on the body zone size Ra in the estimation of the field, the 
multipole moments of the compact bodies, and the equation of motion. Since we introduced 
Ra by hand, it is reasonable to expect that the field does not depend on Ra- We consider 
multipole moments as renormalized moments. It is proved that whatever Ra we take as the 
radius of the body zone, the general equation of motion does not depend on Ra- 

Discarding Ra dependent terms when calculating the field is similar to regularizations 
such as Hadamard partie finie. Then, there is a possibility that at the 3 PN order our 
method might give log^/eR^) type ambiguity that might correspond to the one reported 



in p3|,[13| , |l^| . If such log term should arise we could not discard sRa even practically. 



Our method, however, still differs from other works ||33|jl3| , p^| in some points. Thus, 
besides as a check for the 3 PN equation of motion obtained in |33|JT3|JT^| , it is interesting to 
derive it to study the log type ambiguity. We will be tackling the problem and try to derive 
the 3 PN equation of motion. 
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APPENDIX A: Q¥ AND R\ IJ 



In /ijgn=o; Eqs. ( |3.14|) , and (|3.15|) , the moments J A and appear formally at the 



order e 2l+A and e 2l+2 . Thus Q L A l and R A ^ appear as 

r Lf)Li 

h Bn=0 ~ 1" e 21+1 + 

r A 
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v\lV' J 



h% n =o H f4 f^+i + " 



where we omitted irrelevant terms and numerical coefficients. These equations seem to 
suggest that Q L A % and R L ^ for any I contribute to h^ u at the order e 4 . 

We solve this problem by declaring that h^ u does not depend on the size of the body 
zone boundary, eRa (see Sec. |III C|) . By discarding the terms depending on cRa, we find 



that a finite number of Q L A l and R L A ^ come into play. 

Now, let us consider the n PN order integrand e 2n+2 2n+2^N (where the n PN order 
integrand means the integrand which appears at the n PN order in the evaluation of hg n=Q . 
For example, ^A 1 ^ is 1 PN integrand and qA^ 1 is 2 PN integrand). These integrands, say, for 
the star A, schematically have forms such as 

2n.+2 A ri ( m \ S m2 \ri ( -> -> \ / A 1 \ 

e 2ra+2Ajv ~ e I — I —Aa^a-tiac), (A1J 



m\' m 2 



e 2n+2 2n+2 A^ ~ e 2 " +2 j ^(rfc • n AC ), (A2) 

where the label C denotes the companion star (if A — 1, C — 2). In these integrands, we 
expanded around the center of the star A, and expressed the angular dependence of the 
integrands by the function whose forms are irrelevant here, s — 0, 1, • • - [n — 2] (s = if 
n = 0, 1) and t — 0, 1, ■ • -[n — 1] (£ = if n = 0), (where [a] (a: real number) denotes the 
biggest integer smaller than a). m k and r k are m^rri^ and r^r^ respectively where k, I, and 
n are positive integers and l + n — k. Then from the definition of Q A and R^i Eqs. ( |3.26|) 
and (ET27D , Eqs. (|Al]) and Q give 



/ d \Z-s-1 2ra+2 



e 4 Q^ ~ (efl 



e 4 ^' ~ (ei2 A ) , -*- 1 e 2n+2 , 



where only the e and ei?yi dependence are shown. Now from the argument in Sec. |111 CL we 



are concerned with the terms which do not depend on cRa- Thus when we use the n PN 
order integrand, we evaluate Q S A l3+l1 and R T A H+1% \ 
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To derive h% n=Q at the N PN order, we have to know h^ n=0 up to 0(e 2N+2 ). Then we 
have to use 2n+2^ to evaluate Q A and R A j where n < N. Furthermore, since we can 
discard the cRa dependent terms, we evaluate only Q A and R T A where s — 1, • • -[N — 1] 
and t = 1, • • - [N]. Thus at the Newtonian order, N = 0, Q Lt and R A j do not contribute. At 
the 1 PN order, N = 1, we have to evaluate R A j using 4 [— gt^ L \. Up to the 2.5 PN order, 
we have to evaluate only Q l A , Q\ l , R l i, R\ K \ an d R-a^ using 4 [— gt^ L \ and 6[— 

Explicit calculations required up to the 2.5 PN order are as follows. 

First we calculate 

±- I dS m ^\h^ k y{y{ = AieR^PfC^ + 0{(eR{?+ 2 ), (A3) 

i07T JdB 1 

for 1 even, 

= OdeR^ 1 ), (A4) 



for 1 odd, and where 



Then, Q% (I < 3) becomes 



C L = — i dttn L . 
An 



e *Q% = e 6 £ dS k y L A y\ ( 6 [-gt k L l] - v%[-gt L l}) + 0(e 8 ) 

= le%eR A ) l - 1 P T A 2 v A C^ j + 0{e\eR A ) l+ \ e 8 ), (A5) 



for 1 even, 



= 0(e e (eitU) i+1 ,e 8 ), (A6) 



for 1 odd. Thus Q A does not contribute to up to the 2.5 PN order. The momentum- 
velocity relation becomes 

P l A = P T A v\ + 0^), 

where we omitted the cRa dependence. 
As for R L A ij (I < 2), 
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e 4 i^' 



r ' f dS k y L A y j M [-gt % t L ] + e 6 / dS*3/M ( 6 [-<?4lJ 



+ e 7 / dS ky y A7 [-gt%]+0{i 

JdB A 



e\eR A ) l - l Pl 2 C L k ^ k + 0(e°(eR A ) l -\ e 8 ) 



l-l 8\ 



(A7) 



for 1 even, 



0(e 4 (e^) i+1 ,e 6 (e J R A )' ) e 5 



(A8) 



for 1 odd. Thus R A 3 also does not contribute to up to the 2.5 PN order. 

Finally, since Z A 3 contains J A 3 , while J A J contains Q A 3 , when we evaluate Z A 3 \ we have 
to evaluate not only R A J but also Q A J - From this reason, up to the 2.5 PN order, Q A J and 
Q A ij appear (see Eqs. flOlp and (^32|)). From Eqs. (|AJ) and (P), Qj j = 0(e 4 , e 2 (ei2^)) 
and = 0(e 2 ). These results show that they do not contribute to h^ v up to the 2.5 PN 
order. 

Thus, up to the 2.5 PN order, we have 



> L Bn=Q 



A=l,2 



+ e 



2 M A r A 



+ 0(e b 



fcJUo = 4e 4 E . 

A=l,2 \ Ta 



+ e* 



v A 'r A 



+ 0(e 8 ), 



(A9) 
(AlO) 



where we used Eqs. (|3.21|), (|3.22|), (|3.23|), and (|3.24|) and omitted the eR A dependence. 



APPENDIX B: LANDAU-LIFSHITZ PSEUDO-TENSOR EXPANDED WITH 

EPSILON 



The Landau-Lifshitz pseudo-tensor E51 in terms of h^ u which satisfies the harmonic 
condition is as follows. 

+ \ (sTf* - \g" u g aP ) (g^g* - ^ 7e >% (Bi) 
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We expand the deviation field in a power series of e; 



n=0 



Using this equation, we expand t^ L with e. The point to note is that 1) we raise or lower 
indices with the flat metric rf v and 77^, 2) rf T = — e 2 and r\ TT = — e~ 2 , 3) ^h % \ k = and 
7^ rr ,fc = 0, 4) 5 /T M = 0. Below, vertical strokes denote that indices between the strokes are 
excluded from (anti-)symmetrization. 

I. 0(e 4 ) 



TT 1 

LL\ 



4 [-16irgt 
i[-16ngt? L 
4 [-167r^ L 



0, 

0, 
1 
4 



1 

2' 



(B2) 



2. 0(e 6 ) 



6 [-167T^] = -Z^*^, fc , 

6 [-167T^] = 2 4 h^ M h^ - - A ,h Tk M h^\ 
6 [-167r^ L ] = 4 4 /r li <VV' ? ' ] - 2 4 /i Tr ' (i 4^ )fc ,A ; 



+ l 4 h TT ' (i 6 h^' j) + \ 4 h TT ^\h) k k f - l 4 h TT 4 h TT '\h TT > j 



2 

+ S ij 



411 6 

3 

8 



^4^ ,£6^ — ^4^ 4^ + ^4^ 4^ ,fc4^ 



(B3) 
(B4) 



(B5) 



3. 0(e 7 ) 



r[- 




= 0, 


r[- 




= 0, 


r[- 




1 


~ 2 






1 



(B6) 
(B7) 

(B8) 
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4. 0(e 8 ) 

s[-ifagtY L ] = -\,h^\K\ k 

+ -^4h TT,k 4h l i t k + 4h TT t k4h kl t i 

+ LhT A h TT > k A h T \ k , (B9) 

o 

8 [-16tt^] = 2 4 h T k , H h k[i ' l] + 4h Ti M h kl tl + ^4h Tl M h k k l 

— -4h TT '\h l l tT 

+ 2 6 h T ^\h T \ k - ^ 6 h Tk M h TT >* 
+ 2 G h TT M h T[kA - \ G h TT \h Tk , k 
+ 2 A h TT \h TT M h T[i ' l] + ^ A h TT A h TT '\h T \ 

- \ A K k A K\ kA h^ + ^h Ti 4 h TT ,i4h TT \ (BIO) 

s[-l6n 9 q L ] = l - A h l \h TT ' k A h TT k 

- 2^\h^\ x + \ji TT ' { \ti\f - 2 A h TT > { -\h? ) \ l + ^ TT ' (i 6/i'V j) 
+ 2 6 h Tl ' (i A h j) \ l - 2 6 h T(tM 4 h 3) \i + 2 A h Tl > { \h j) \i - 2 6 h T l ' (i A h^' j) 

- 2 A h Tl '\h\ T 

+ 4h lk ' l 4hh,i + 4h ik M h jl tl - 2 A h kl ^\h?\^ 
+ l -4h kl \h kl ^ - - A 4h k k \h\i 

+ i 4 /i Tr ' (l 8^ |Tr|j) 

+ A h TT A h Ti M h Tj ' 1 + A h TT A h rl \h T ^ 

+ 2 A h T1 A h TT ' (i A h j)k >k - 2 A h TT A h rk > { \h j)T , k 

+ A h T kA h Tk ^ A h lTTl ' j) + ^ A h Tk M h TT ^ A h j)T 

- l - A h T \ kA h r ^\y )k - l - A h rT A h k k >\h}^ 
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h T \h TT ' ( \h^' j) - - 6 h TT 4 h TT '\h TT ' j 



— ±lb 4 

+ -Lh TT )\h TT '\h TT ' ] 
4 



+ 5 lJ ~4h TT \h T ^li T ,/ — 4 /i TT \h TT \ k Ji j — 4 /i TT \h T ' ih T [k fy 

_o 

i i 

1 -7-1- 1 T-T- 111* 

r - 2' 



+ - Ah TT Ah TT)l Ah k ki — -Ah Tk Ah TT ' 1 ih T ki 



— - 4h Tk \h TT t kAh Tl ,/ + -4^ H 4^ rT ,fc4^ TT ,/ 

4 o 
+ 4,h T k,uh M ' T — -4^ rfc ,fc4^«,T 

1 1 1 

+ -4fl ±n-km,l — ^4^ 4^fcZ,m + g4^ fc.W'i m 

3 k t kl 

— ^h T )k §h T j + 2 A h T ' Qh T [kt i] 



+ §h TT \uh lk t k — ^6^ rr,fc 4^'i,fc 

+ -,h T \h TT M h TT ' k + Lhr^hr^hT'" 
4 2 

- luh TT )\h^ k 4 h T \ k 



— - 4 Al 8 ft >fe 



1 

8 f 



~6 ft 6 ft ,fc 



(Bll) 



5. 0(e 9 ) 



9 [-16tt^] =0, 



9L -167n?*li] = 2 4 /i rr ,, 7 /r [M] - 74^5/^ + ^h TT '\h TT 



1 

4*" ' J '""' ' 4' 
[-16tt^] = ^ Tr ' ( ^ M ' j) + l^'Vi^ _ l^hT^hr^ 

+ l,h^ k ,h T \ k5 h^ - \ 4 h TT M h TT ^^ k 

o z 

- 2 4 h T ^\h j \, T 

+ 2 A h TT ' (i 7 h j) \ T 

- 2 4 /iV ( > |rfe| " ?) + 2 4 /i r(l | fc |7^ |rfc|,i) 

- 2 4 /i T ( i | fc | 7 /i J ' )T ' fc + 2 4 K r W 7 hP>, k 



(B12) 
(B13) 
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+ 5 tJ 2 4 /i T ■ 7h T [ k j] — 4 /i TT \k7h T )T + Ji T k,ibh .t 
1-3 l kl 

4 4 o 

1 11" 

I UTT,k ITT ITT UTT,k ll W~T,k ITT 

+ ^4^ 4^ ,fc7 ft _ ^4^ 7 ft Z,fc ~ ^4^ 9 ft ,fc 



(B14) 



APPENDIX C: SPIN-ORBIT COUPLING FORCE 

It is well known that definition of a dipole moment of the star, which we equate to zero 
to determine the center of the mass of the star, affects form of the spin-orbit coupling force 



e.g., @])- I n the paper I we have chosen d\ = 0, where 



D\ = d A + e 2 M^v\ + 0^). 



Instead, in this paper we choose D A = 0. The corresponding spin-orbit coupling force is 
( dv\\ A 2V k 



mi 



dr 



—e 



so 



12 



m x M$ + m 2 hlf) A lk + {m x M lk + m 2 M[ k ) A* 



+ e- 



where A. lj = 5 lJ — 3n\ 2 n{ 2 . 



rf 2 



m 2 M lk v\ - m x M l 2 k v k 2 



lk„.k 



A' 



APPENDIX D: DEFINITION OF THE MASS 

In Sec. |jVl we defined the mass of the star A as 

m A = \\mP T A . (Dl) 

We shall explain the meaning of this definition. Also we shall explain the reason we define 
the mass as above. 

Let us consider the origin of the e dependence of P A . The origin can be classified into 
three categories: 1) The effect of the motion of the star A itself and the gravity of the 
companion star. 2) the eR A dependence which comes from our splitting of space time, and 
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3) the deviation of the scaling law which we impose on on the initial hypersurface. We 
shall explain them respectively. 

1) This type of e dependence of P\ can be specified by solving the evolution equation 
for P A . In fact, for example we can rewrite Eq. ( |4.5| ) as 

pr . 2 1 2 . 2 m i m 2 . 2 2m 1 m 2 4 

P{ = m 1 + e -m 1 v 1 + e h e h C(e ). 

2 r 12 r i2 

This is the Newtonian energy (times \f—g) of the star 1, if one lowers the upper index r 
of P[ down with Newtonian metric and multiply it minus sign. At least up to the 2.5 PN 
order, this interpretation really works (see the appendix [H]). Thus we can specify the type 
1 e dependence of P A by solving the evolution equation for P T A functionally. 

In fact, our definition of the mass, Eq. ( pi|) , comes from this observation; if the star A 
stayed at rest and there were no companion star, P T A would be the ADM mass of the star A. 

2) Since we defined P T A as the volume integral over the body zone, it must depend on the 
size of the body zone, €Ra- As we have stated in Sec. [Ill C| , we ignore this €Ra dependence. 
But even not doing so, when we let e go to zero, this eRa dependence disappears. Let us 
explain how this occurs. 

First, because we assume a non-singular source for the star, the gravitational field on 
the body zone boundary must be smooth. In particular we can use the post-Newtonian 
expanded gravitational field near and just inside the body zone boundary. Thus we can 
estimate the cRa dependence using the post-Newtonian expanded gravitational field. Then, 
the most dangerous term (,i.e., the term which seems to diverge when e goes to zero) in the 
integrand of the volume integral Eq. ( |3.37|) at the n PN order has a form (see the appendices 
0and0) 



Alf ~e 2 " +4 /rt,ft ,~ 



n-l 



m A m A 



where we omitted the indices. Then we can estimate the cRa dependence at the n PN order 



as 

1 A) \eR A ) {eR A Y 
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An explicit calculation (see paper I) shows that up to the 1 PN order (,i.e., n < 1 in the 
above equation) the cRa dependence takes a form like e 2 m A /eRA- Thus the cRa dependence 
vanishes when e goes to zero. 

3)Though we impose scaling law on T^ u as Eqs. (|2.2|), (|2T3|), and (|2.4|), there must be 
deviation from these scalings in the integrand A^ . 

One reason of this deviation is that since the velocity of the spinning motion of the star is 
set to be 0(e) (slow rotation), the scaling law on the initial hypersurface can not be exactly 
satisfied by A^f . Besides this, the possibility of obtaining a sequence of solutions which has 
the scaling imposed in this paper has not been established. We only assume the existence 
of such a sequence. 

The deviation may also come from the evolution effect. For instance, as the separation 
of the binary gets shorter and shorter, tidal effect stimulates oscillation of the constituent 
stars. Such an oscillation may let the stars not satisfy the scaling law. 

If we can construct a sequence of solutions which has the scaling law (Eqs. (|2T2|), ([2.3|). 
and Q2.4p ) at least as their leading order, and when we consider a binary system whose sepa- 
ration is long enough, then we can incorporate the above three deviation into our formalism 
by expanding the density pa in an e series (Let us use the pa instead of A^ for notational 
simplicity) . 

Pa = e~ 4 (-4)p A + e _3 (_3)PA H • 

From this expansion we have an expansion of the mass in the e series. 

m A = (o)m A + e m m A H , 

where 



(n) m A = lim^ d 3 y (n _ 4) p A - 



>b a 

Here we need to take e zero limit to describe the star A as a point-like particle. 

Thus the e zero limit in Eq. ( pi[ ) is taken to 1) remove the effect of the motion of the 
star A itself and the effect of the existence of the companion star, 2) remove the dependence 
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of rriA on the body zone boundary, cRa, and 3) achieve the point particle limit, while not 
removing the effect of the deviation of from the scaling law imposed initially. 



APPENDIX E: THE 2.5 PN GRAVITATIONAL FIELD 

We list the deviation filed W up to the 2.5 PN order. We have checked the harmonic 
condition up to the 2.5 PN order; <jh' IU )V = 0(e s ). We can transform h^ v into the metric by 
the formulas given in the appendix |G]. The resulting metric perfectly agrees with the result 
in 



12 



Here are some notations. 



f (lnS) = i(~ r i 2 + 3r i r i2 + r u - 3rir 2 + 3r 12 r 2 - r\) 
3o 

+ -^{rl-r 2 12 + rl)\nS, 

= T^{- r l ~ 3rir 12 - r\ 2 + 3r x r 2 + 3r 12 r 2 + rf) 

lo 

+ \{rl + r 2 12 -rl)\xiS. 



(El) 



(E2) 



They satisfy A/( ln5 ) = lnS 1 and A/^ 1, l > = ri/r 2 . These two potentials (and other useful 
inverse Laplacian formulas) can be found in the appendix of JTTJ . 



pr 

4e 4 E — 

A=l,2 T A 



a,b=l,2 



ran ru A %2 r A ' dr 2 u=1;) 



4 7 d 
3 dr 



E Plvl - 



, A=l,2 



pr pr 
T\1 



■ M E 

r\// 



Pt3 dt2 



A=l,2 r -4 .4=1,2 r A 



E -^(s(n A -VA) 2 -7v 



dr 2 



2r 12 



S 
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+ 6 S 



4 d \m 1 m 2 ^ / 2 _ n2 \ 5m 1 m 2 (m 1 + m 2 ) 

" 3( " 12 ■ ) + 5? 

5mim 2 _ _ 13mim 2 1 
2r 12 ' Ul )( ni2 ' ^ 2^ 2 _ ^ ' ^ J 



.2 
12 



0(e iU ). (E3) 



pr2 pr 
-^1 -*2 



8 3 173 

3 3rir 2 ; 



+ 



9ri 



+ 



r ,3 
12 

80 



52 1 

~0 Q 



103 



23r 2 



+ 



V3riri 2 r x r 2 
-7 





r?2 


O r 2 „3 
z ' 12' 2 


2r 2 r 


15r 2 


8r 2 


2r?r 2 2 


r 3 ri2 


2 


16 \ 



+ 



21 



4r 2 



31 2 

3r i2 r 2 riS" r^S r 2 S J 



+ 



1 v rir 2 S 2 riS 
7 2 2 



(ni • v ± ) 2 - 



12(ni -ui)(ni2 -ui) 



^12^2 



S 2 r l2 S 



(n l2 ■ vif 



S 2 



( 28 

+ ( 



+ 



28 



S 2 

16 2 

+ 



+ 



16 



+ 



3rir i2 r\r 2 3r 12 r 2 nS r 12 S r 2 S 
24(n 1 .u 2 )(n 2 -tf 1 ) /32 32 \ ^ ^ 



S 2 




-16 


16 \ 


S 2 


r 2 Sj 


) (ft ■ 


v 2 ) + 



8r2 



2rf 2 r 2 ' 3r 2 2 r 2 2r\r X2 r 2 2r\r\ 2 ' 2r\r\ 2 rfr 12 r\r\ 2 r\r\ 2 ' rfrf 2 ) 



(n 2 ■ vif 

16(ni • v 2 ){ni 2 ■ vi) 
S 2 



S 2 
2 



10gi ■ vi)(n 2 ■ v 2 )' 

s 2 



12(gi ■ t7i)(wi2 • ff 2 ) 
S 2 

+ (1 <-> 2) 



(E4) 



PaVa 



A=l,2 r ^ 



+ 7<} + ^p-( W f + v k 2 ) (S ki - n{ 2 n k 2 ) 



A=l,2 



Sr 12 



+ 8P[P 2 T (v{ + vi 

JDT pi 
r l r 2 



1 



r 12At^2 r ^ r l r 2 



- 16 
+ 12 



s - 2 K« - n[)(n k 12 + n k ) + v k (n k 2 - n k )(n\ 2 + n 2 )} 
{vi(ni2 ~ n1)(n[ 2 + n 2 ) + v k (n[ 2 - n[)(n k 12 + n k )} 
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d 



2 E p( p >w) 



A=l,2 



(E5) 



fcii = 4e 4 ^ 



A=l,2 



ta 



-2e 



A=l,2 

d 3 



pr2 



-n 12 n 12 



8P 1 r P 2 T 
S 2 



(^fc^'i - (nia - ni) (fc (ni 2 + n 2 )° 



dr 3 



+e6 E V 

A=l,2 



E Pa z A z A 
A=l,2 

API 2 



5 lJ 



v\ + (n A ■ vaY 



v 2 A n\n A - ^— ' ^ 2n A n A + (nA ; ^)nVj + jj^>j ) 



+e 6 6 if ij ' 



+2e e 



<9r 2 



E ^>>>a + - E *W - 3^>i) 

A=l,2 60 A=l,2 



A=l,2 



S*3 

, + —PIPo In 
K ej 3 1 2 



5 



i ( - -<f^ ) p; p. 



2U/e 
T fd 2 f^ 5 kl 



dz{ k dz l 2 ] 6 



In 5 



7 <9 3 

+0(6 8 ). 



-q E m AV l A v A r\ + 



A=l,2 



7mim 2 
ri2 



22 



E z aA - ^m x m 2 r X2 b i: > + T ^ 1 n\ 2 n\ 2 ^ 



A=l,2 



A=l,2 



(E6) 



Here g^P- 7 is 



, r2r)T , , 23 3ri 163 11 17 r 2 17r, 



Ap T2 p T J __ : 1 _j 1 I 

- ! ' " ' 2Ar\S Ar\ 2 S 2Ar 2 2 S l2nr 12 S 8nr 2 S 8r\ 2 r 2 S 8r 2 2 r 2 S 8r 12 r 2 S 

llrg 17r 2 23r 2 llrj _ 23r 2 23r| \ :J 



12r 3 2 ,S 8r!r 2 2 5 24r 2 r 12 ,S 12r ir f 2 S 2Ar\r\ 2 S 2Ar\r\ 2 S , 

+ ^ i_L_ _rf_\ 4 ■ /^l 26 r,_ '_ _J_\ (i ■) 

+ l v 24r 3 2 24r 2 r 12 + 24r 2 r?J 11 + ^3nr 2 2 3r 12 S 2 3r 2 2 S 2 r 2 2 S j 1 12 

/ n 9 2r 2 15 

+ 77 o 7777 H 7777 + 77-^777 + ^r^7 I ni,n 



,3r 2 2 S 2 r 12 S 2 3r 2 2 S 2 r 2 2 S 
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+ [ ~ r2 + Vl + ri - 17 - r 2 ] n ^J) 
\3r 2 12 S 2 Qrf 2 S 2 3r 2 2 S 2 2r 12 S 2 6rf 2 S 2 I 1 2 



-n 25 \ 



+ 4P^ j-2, 2 - \5n-) + 8 - |^)< 

_ gxUyy d 2 lng , H <9 2 ln,S 
, <9 2 ln,g 3 <9 2 ln,g \ , , 



1 \ -94^(1-1) 



2" " J dzi k dz l ldz?dz? 1 1 J 
+(l<->2). (E7) 

Here in the terms (1 <-> 2), f^ 1 ' 1 ' is replaced by p 1 ' 1 ', and this function satisfies A/^ -1 ' 1 ) = 
r 2 jr\. Its explicit form is the same as Eq. (E2) but with the labels 1 and 2 exchanged. 

APPENDIX F: X PART 

Since Xtv^ ',a/3 itself is conserved, i.e., XaT Q/3 ' ,aPn = 0, the P£ by itself is conserved, here 
we define Pj£ as 

n x = e-'f Axrv (fi) 

Jb a 

First we derive the functional dependence of P^ on m^,^, and ri2- By the definition 

of 

l^X TTaP , a /3 = {h Tk h Tl - /T^),^ 

167Tx" a/3 ,a/3 = (/i TT ^ fc - h Ti h Tk ), Tk + {h Tk h U - h Ti h kl ) 



Mi 



thus, we can obtain the functional expressions of P£ using Gauss' law. To obtain these 
expressions up to the 2.5 PN order, we need Ji^ and 5 ft/ 4 * '. See the Sec. |V|. The results are 
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pT _ A m x m 2 
lx ~ 3r 12 



4y2 + mi _ 2mi 
na r i2 



- e 5- mi (3)/^ + o(e 6 

3 A x = ^a + o(6 6 ). 



(F2) 
(F3) 



Now we move on to the equation of motion for P% x - 



d P A x 



—e 



dr JdB A 
Evaluation of surface integrals gives 



f dSkxff^ , a f3 + e i v\i dS, 

JdB A JdB A 



kX N ,a/3- 



(F4) 



dR 



T x _ 4 mim 2 (ni2 • V) 



dr 



dP l 
dr 



—e 



3r 2 



12 



4V 2 + — m2 ' 4mi 



4 mim 2 (ni2 ■ V) j 
£ z o f i 



3r 2 



12 



+ 
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+ 

ri2 ri2 
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3rf 2 



3 £ 



12 



4y2 + _ 2mi 

^12 7-ia 

i(4)W ^^12(3), 
"H^l ^orbi 2 J < 

r 12 



(F5) 
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These two equations, Eqs. ([F5|) and (|F6|) , are obviously consistent with Eqs. (|F2|) and 



(ED- 



APPENDIX G: METRIC EXPANDED WITH EPSILON 



Here we give the metric expanded in an e series for convenience. 





1 + 


-e 4/1 


+ 




6 /l TT — Jl k k — -(4^ TT ) 2 


+ 
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0(e 6 


) 



>h TT + 4 h k k --( 4 h TT ) 2 



+ ^-[ 7 h rT + 5 h k k ] 



(Gl) 
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+4 
+4 



3 1 5 

ih TT + §h k — -eh TT 4: h TT — - A h kih TT + Ji T ^h 1 \ + -U^ TT ) 



1 



)h TT + 7 h\ - - 4 h TT 7 h TT - -ih TT 5 h k 



+ 0(e 6 ) 
e 9ri — —e 411 — e 
e~ l g tl = S ij 



e 5 7 h Ti + 0(e 6 ) 



(G2) 
(G3) 



1 + e 2 ^ 4 h TT 



+ t 



1 

2' 
<^ 
2~ 



+ 0(6 6 ) 



(G4) 



APPENDIX H: MEANING OF P T Ae 



In this section we explain a meaning of P A q- First of all, we expand in an e series the 
four velocity of the star A normalized as g^ u u A u A = — e~ 2 , where u\ = u T A v A . By the help 
of the appendix [TJ, we have 



+ e 



1 + e 

1 

4 



e 5 - 
4 



4 



5'i fe 
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0(6 6 ). 



(HI) 



This is a formal series since the metric derived using the point particle description diverges 
at the point A. Now let us regularize this equation by letting zero while simply discarding 
the divergent terms. For example, by this procedure <eh TT becomes (see Eq. ([4.6|)) 

I -To 



[<6/r r r* 



4e 



4 J 2 

ri2 



m 2r/ ^ . 2 2 mim2 



-2^{(ni 2 -^r-^}-16- 



+ 7^ 



(H2) 



r 12 '12 '12. 

for the star 1. In the above equation, [f} e A xt means that we regularize the quantity f at the 
star A by simply discarding the divergent terms 
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Evaluating Eq. ([HI]) and (|G1|) by this procedure, then comparing with Eq. ( |7.2|) with 
Eqs. ( |6 . 1 6| ) and ( |6.17| ), we find at least up to the 2.5 PN order 



P T A e=m A [V^u\]T. (H3) 

It is important that even after discarding all the divergent terms, there remains non-linear 
effect. See Eq. (1617)) and Eq. 



Eq. ( |H3| ) is natural. And note that we have never assumed this relation in advance. 
This relation has been derived by solving the evolution equation for P\ e functionally. 

The result of the above procedure, "discarding all the divergent terms", coincides with 
the result of some regularization scheme, such as Hadamard partie finie, at least up to the 
2.5 PN order. We stress, however, that this coincidence is true for the calculation of P\ 
up to the 2.5 PN order. This happens because only 4/1*** ' , ^h l \ fh TT , and §h TT are needed for 
the calculation of the 2.5 PN P\ (see Eq. (|H1|)). When one calculates the 3 PN P T A , for 
example, %ti rr is needed. Then "discarding all the divergent terms" and Hadamard partie 
finie give different results (to see this, regularize H\ by both ways, for example). 
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